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Abstract 

Since Davidson patented Circulation Control Airfoils in 1960, there have been 
only two aircraft designed and flown with CC. Designing with CC is complex for the 
following reasons: the relation between lift increase and blowing momentum is 

non-linear; for good cruise performance one must change the wing geometry in 
flight from a round to a sharp trailing edge. The bleed air from the propulsion 
engines or an auxiliary compressor, must be used efficiently. In designing with CC, 
the propulsion and control aspects are just as important as aerodynamics. In this 
paper these design aspects have been examined and linearized equations are 
presented in order to facilitate a preliminary analysis of the performance potential 
of CC. The thrust and lift requirements for take-off make the calculated runway 
lenght very sensitive to the bleed air ratio. Thrust vectoring improves performance 
and can off-set nose down pitching moments. The choice of blowing jet to free 
stream velocity ratio determines the efficiency of applying bleed air power. 


Introduction 

Davidson (1960) patented the Circulation Control concept. The initial application 
was for cylindrical airfoils. Kind et al (1968) provided experimental data for the 
elliptical airfoil. The ellipse is desirable for helicopter blade applications because it 
has some lift generating ability in case of power failure and is structurally rigid. In 
recent years, it's leading and trailing edges have been modified and camber was 
added to improve the lift characteristics. More wind tunnel data are available on 
the basic and modified elliptical shape CC airfoil than all other configurations 
combined. The elliptical shape with it’s maximum thickness near mid-chord has a 

center of pressure near mid-chord and thus a steeper adverse pressure gradient 
than conventional airfoils. This thickens the boundary layer upstream of the 
blowing slot and renders the Coanda turning efficiency very sensitive to geometry, 
angle of attack, Reynolds number and turbulence level. Various investigators have 
found significant differences in the lift-to-blowing momentum augmentation ratio for 
seemingly similar configurations. The non-linearity of this ratio with the blowing 
coefficient further adds to the complexity of selecting an optimum CC configuration. 
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When a CC rounded trailing edge is added to a conventional high lift airfoil, 
complete with a leading edge modification to prevent separation, it’s performance is 
much more predictable. Such airfoils have their maximum thickness further forward 
than the ellipse and therefore possess a less steep adverse pressure gradient than 
the ellipse. The reduced pressure gradient, together with the usually greater 
chord length and Reynolds number, minimizes the difference in boundary layer 
shape ahead of the blowing slot for CC modified CTOL wings. Such airfoils can be 
used for STOL aircraft and provide more comparable lift augmentation ratios and 
linearity with angle of attack. The greater thickness of a STOL aircraft wing as 
compared to a helicopter blade, permits incorporating advanced leading edge devices 
such as Kruger flaps, slats, or drooped leading edges with or without blowing. 
L^rge ducts fit inside a STOL wing which permits the power efficient use of blowing 
au- at medium pressure and temperature. An ejector built around the blowing slot 
can provide boundary layer suction just upstream of the rounded Coanda jet 
turning surface. For added lift, the chord length may be increased in the CC mode, 
however the' rounded trailing edge must be retracted in order to obtain efficient 
cruise performance with a sharp trailing edge. The availability of flaps in 
conjunction with CC is desirable for drag control during descent, adding flexibility 
and safety in the case of a blowing air failure. Even though the STOL CC airfoil 
geometry is much more complex than the modified ellipse, it's aerodynamic behavior 
and linearity is similar to that of the conventional CTOL airfoil. This is evident 
from the wind tunnel model test data by Englar (1975) on various STOL airfoil 
configurations. The flight tests on the only two CC STOL aircraft built and tested 
showed good agreement in lift augmentation ratio and performance, even though 
these aircraft were significantly different in geometry and wing loading. 

Since the CC air mass flow requirements are low they may be provided by a 
relatively light weight auxiliary turbo-compressor or a compressor driven by shaft 
power take-off from the thrust engines. In this manner the blowing pressure is 
independent of the thrust level, thereby simplifying the design and operation with 
CC. However, the weight, cost, and reliability penalty of an additional compressor 
makes the use of jet engine bleed air more attractive. The thrust loss associated 
with compressor bleed air is a complex function of bleed pressure and throttle 
setting. The bleed mass ratio is one of the most important parameters in 

determining the effectiveness of CC during the take-off and landing phase. 
Minimizing the blowing air ducting heat and pressure losses, as well as the 
throttling loss across the bleed air flow control valve, greatly effects the ratio of 
blowing thrust generated to engine thrust lost. With CC, the blowing coefficient is 
in general less than 5 % of the lift coefficient, therefore the blowing momentum is a 
small fraction of the aircraft weight. Most of the momentum is used to energize the 
boundary layer along the Coanda turning surface. The small fraction remaining in 
the wall jet after separation is often negligible relative to the take-off thrust 
coefficient. Even in the high drag landing configuration, the effect of Cp on the 
reduction in the drag coefficient can often be ignored. 


WVU CC STOL Demonstrator Design Aspects 

In 1968, the Office of Naval Research awarded a research contract to West 
Virginia University to investigate the theoretical and experimental aspects of 
Circulation Control. This was about the same time that Williams (1970) at NSRDC, 
started to investigate the performance of elliptical airfoils to be used on a high lift 
helicopter. As a result, significant improvements were made in the lift-to-blowing 
momentum augmentation ratio. However, it became apparent that this ratio was very 
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dependent on the model aspect ratio, slot adjustment mechanism, jet to free stream 
velocity ratio, tunnel blockage factor, turbulence level and Reynolds number. To 
eliminate tunnel wall interference and achieve high Reynolds numbers, WVU 
researchers decided to flight test a CC wing on a light aircraft. The first CC 2-D 
STOL airfoil built at WVU had a 4 ft chord and incorporated a drooped and blown 
leading edge. A low pressure blowing slot produced an air jet tangentially to the 
cylmdricai flap hinge. The airfoil could be operated in the blown flap mode or as a 
CC airfoil after the 15% flap was folded fully forward and flush with the underside 
of the airfoil. In 1970, this WVU model was tested in the 8 x 10 ft NSRDC tunnel. 
^ en though it performed as anticipated, the following improvements were made. 
The flap was folded out in the CC mode to increase the chord instead of decreasing 
it by 15%. The rounded Coanda surface remained stowable in flight, to provide a 
sharp trailing edge for low drag cruise. A choked flow nozzle was selected for 
better span-wise blowing uniformity and slot adjustment screws were eliminated. A 
built, in ejector was designed to improve the jet thickness, temperature and Coanda 
turning of the otherwise supersonic jet. The ejector entrained air provides needed 
structural cooling and boundary layer suction at the flap hinge. Reducing the 
boundary layer thickness upstream of the blowing slot greatly improves the Coanda 
turning. The drooped leading edge was designed by Norio Inumaru, program 
manager of the Japanese USB QSTOL, and proved to be so effective, that the leading 
edge blowing could be eliminated. In 1971, the improved CC STOL airfoil (figs. 1 
and 2) was tested in the WVU wind tunnel, Loth (1973). At the same time 
construction was started on the WVU CC STOL Demonstrator Aircraft, using materials 
from a BEDE 4 home builders kit. A 200 HP GTC 85-72 APU was selected to provide 
up to 1 kg/sec of air at maximum 2 atm gage, to the CC blowing slot. The blown 
ailerons could be drooped. Inboard and outboard fences were installed for both 
aerodynamic and structural reasons (fig. 3). Both chord-wise and span-wise 
pressure taps were installed on the wing. A sliding blowing air dump valve was 

installed to provide direct lift control. The DLC valve was actuated by a lever on 
the throttle quadrant. A splitter valve, linked to the aileron controls, provided 
optional roll control by differential blowing (fig. 4.) 

On April 10, 1974, pilot Shawn Roberts (1974) started 25 hours of flight testing 
to determine the performance potential and handling qualities of CC. The WVU CC 
STOL Research program was summarized by Loth and Fanucci (1974). The lift 

coefficient appeared to increase linearly with angle of attack and with the square 
root of the blowing coefficient. At high blowing rates, the trim lift coefficient was 
considerably lower than it was on the CC flap. The loss in lift is due to the 
ownload of the stabilator, associated with the large nose down pitching moment. At 
c Lmax> with thrust power at idle, the CC blowing at C M = 0.12 increased the flap lift 
coefficient from 2.1 to 5.2, or with a lift augmentation ratio: ^C^/Cp z 3.1/0.12 = 25.8. 
At the same time the trim CLmax increased only from 1.98 to 3.8, or with a lift 
augmentation ratio: *C L /Cp = 1.82/0.12 = 15.2 and C M = 3% of C Lmax . This is 

considerably lower than the 2-D augmentation ratios obtainable with an elliptical 
effectiveness is better described by the proportionality constant Cr r 
a Cl/(C p) u,D , which does not change with Cp for a given airfoil. 

Englar (1978) published wind tunnel data on a three dimensional CC model of a 
modified A-6A Navy Cruisader. Even though the CC wing configuration was entirely 
different from the WVU model, the lift augmentation ratio's obtained agree very well. 
His date shown in figure 5 are nearly linear with angle of attack and with the 
square root of the blowing coefficient. The solid lines represent an empirical curve 
fit with dC L /d« = 4.74 and C B = 6.1. In 1979, Englar’s efforts resulted in 10 hours 
of flight testing of a CC modified A-6A Navy Crusader. All the flight test data, Carr 
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(1986) and wind tunnel model data, appear to have about the same constant of 
proportionality (Cb), as shown in Table 1. 


TABLE 1 COMPARISON OF 3-D STOL WING TEST DATA 


SOURCE OF TEST DATA 

| BLOWING) 

TAIL-OFF, WING ONLY 

TRIhWED 

AIRCRAFT) 


1 C M 1 

*Cl | 

CB = 

*c L 

|CB 

=f£i | 


1 1 

1 

y/cJT 

1 

1 

■JCfT * 

1974 WVU CC STOL Aircraft 

1 1 

1 



1 

1 

Flight Tests, Loth (1974) 

1 1 

1 



1 

1 

a) Based on cruise wing area 

1 .12 1 

2.3 1 

6.6 

1 1.8 

1 

5.1 1 

b) Based on CC wing area 

1 .10 1 

2.0 1 

6.1 

1.5 

I 

4.7 1 

1979 Grumman A-6A CC STOL 

1 0.05 1 

1 


1 1.3 

1 

5.8 1 

Aircraft Flight Tests, 

1 0.025 1 

1 


1 0.85 

1 

5.4 1 

Carr (1986) 

1 1 

1 



1 

1 

Wind Tunnel A-6A Model Tests 

1 0.20 1 

2.75 1 

6.1 


1 

l 

Tail-Off, Englar (1978) 

1 0.10 1 

2.10 1 

6.6 


1 

1 


1 0.05 1 

1.36 1 

6.1 


1 

I 


1 0.025 1 

1 1 

0.90 1 

1 

5.7 


1 

1 

1 

1 


The influence of CC blowing on the lift curve slope (dCi,/da) is negligible. 

The effect of increasing the wing chord is the CC blowing mode by folding 
out a flap, is to multiply the lift coefficient by the ratio of the wing area 
increase. When a flap is used in conjunction with CC blowing, the lift 
increases as a function of the flap deflection angle (<5-f). This may be 
estimated from thin airfoil theory as: 

dC L /d« f = 2n - 2e f + 2 sin (e f ) (1) 

where the flap hinge locater angle (flf) is: cos 1 ( 1 — 2 cf/c). 

Combining all the factors contributing to the lift coefficient with 
CC: 


_ c (with CC ) 
C L c (without) 


a 



+ 



+ Cg J Cm 


( 2 ) 


Linearized CC Lift Augmentation Ratio. 

With a choked flow isentropic CC nozzle, the (2-D) blowing momentum is given by: 
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m .V . 

J J 


P .hV 2 
J J 
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m.V. 

C 

p q c 



V .1 2 

v“ 



(3) 


For a specified blowing pressure ratio and slot height (h) the momentum is not a 
function of the temperature (tj). The reason is that the exit density is proportional 
to (1/tj), whereas the square of the velocity is proportional to (tj). In the case of 
regulated and fixed bleed air flow rate, any heat loss in the duct lowers the nozzle 
total pressure and therefore the blowing velocity (Vj). Assuming an isentropic 
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compressor and duct, the jet exit static temperature and density would be identical 
to the ambient values. In practice, the increase in temperature (tj) due to 
compressor inefficiency is approximately off-set by the duct heat and pressure loss, 
thus the density ratio in the blowing coefficient may often be ignored. 

In Figure 6, are wind tunnel test data by Englar (1981), plotted as a function of 
the lift coefficient versus blowing pressure ratio. The result is highly non-linear. 
However, when the same data are replotted as a function of velocity ratio (fig. 7), 
they fall on a straight line! The data appears to fit a single empirical equation with 
the exception of the data for a very narrow slot operating at high velocity and 
pressure. The reason may be that the slot deflection under high pressure is more 
significant at small slot heights. Another explanation may be that at low mass flow 
rates the expansion becomes more isothermal which would tend to increase the jet 
velocity, when testing with room temperature compressed air. It is noteworthy that 
velocity ratios below 1.0 seem to provide no lift augmentation, thus there is a 
minimum value (C^ = 2h/c) below which there is no lift augmentation. Such low 
values, typically below 0.002, are never used and this term becomes negligible at Cp 
greater than 0.02. This allows one to linearize the lift augmentation ratio with 
respect to the square root of the blowing coefficient. In figure 6, the coefficient of 
proportionality (C B ) equals: 9 for h/c = 0.0003 and 10.9 for h/c = 0.0012. Note that 
the magnitude of (Cg) for a three dimensional wing is considerably lower and closer 
to 6 for the conditions shown in Table 1. 

Because the engine thrust is reduced when blowing power is extracted, the 
choice of the blowing velocity to flight velocity ratio, becomes very important. At 
constant (C^), the required blowing power (Pb) decreases in proportion to the 
reduction in velocity ratio. Then for the linear lift augmentation model, the lift 
coefficient ratio ( A Cb/Pb) increases with decreasing velocity ratio, as shown in 
figure 8. However as (Vj/V) approaches 1.0 the value of aCl reduces to zero, this 
makes the constant power curve peak at: Vj/V = 2.08. A similar optimum velocity 
ratio was found with boundary layer energization by tangential blowing to prevent 
separation in an adverse pressure gradient, Boasson (1985). Note this velocity ratio 
optimization does not consider blowing air duct characteristics such as: losses, size 
and weight. When these factors are incorporated the optimum velocity ratio is 
usually greater than 3. 

The effect of duct losses in heat and pressure, are shown in figure 9 in terms 
of percent blowing momentum recovered at the nozzle. This example is for a typical 
case where the compressor bleed air is extracted at: 11 atm and 375* C. As can be 
seen, a 50% loss in air temperature is more detrimental than a 50% duct friction 
pressure loss. The loss in temperature will be noted by an increase in pressure 
drop across the valve controlling the bleed air mass flow rate. The reason is that 
at constant mass flow rate, the loss in temperature results in a reduction of the 
nozzle total pressure. The reduced nozzle pressure is then transmitted through the 
subsonic ducting to the control valve. 


Thrust Loss Due to CC Blowing Power Extraction 

The air power requirements for CC, by blowing over a rounded trailing edge, 
are greater than for conventional boundary layer control, however less than for 
any other type of powered lift system such as: the jet flap, augmenter wing, USB, 
etc. The CC 3-D blowing momentum coefficient, based on the entire wing area rarely 
exceeds 5% of the lift coefficient, or at Cl = 6 find C^ < 0*3. Thus even in the 
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absence of thrust vectoring, the blowing momentum is less than 5% of the weight. 
For most STOL aircraft this means less than 10% of the thrust. The optimum CC 
blowing velocity is at least three times higher than the lift-off velocity on take-off 
and of the same order of magnitude as the cruise velocity. The propulsion system 
exhaust velocity for a propeller driven aircraft is also of the same order of 
magnitude as the cruise velocity. Therefore, the jet-kinetic power at the CC nozzle 
is usually less than 10% of the available thrust power. For jet engine propulsion 
the exhaust velocity is at least double the cruise speed and the CC jet-kinetic 
power at the nozzle is usually less than 5% of the thrust power. 

An auxiliary compressor for CC blowing can be driven efficiently with a power 
take-off from propulsion engines driving a propeller. Propeller aircraft suffer a 
reduction in thrust with forward speed. The power extraction for CC blowing is so 
low, that the anticipated thrust loss is roughly offset by the associated reduction 
in lift-off speed, so that lift-off with, or without CC, occurs at about the same 
thrust to weight ratio! The relative thrust loss due to power take-off can be 

computed using actuator disk theory and a coefficient of performance (C p ). Without 
power take-off the static thrust is defined by: (T 0 ), the induced velocity at the disk 
by (w Q ) and in the wake by (2 w Q ). From the change in momentum of the mass flow 
rate (m) through disk area (Ad) find: 

T = m 2 w =2 pA.w 3 (4) 

o o do 


The ideal power is related to shaft power (P s ) by 

d o 
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(5) 

(6a) 


T ~ 
o 



(6b) 


With forward velocity (V), the induced velocity reduces to (w) and the thrust to (T). 
Assuming constant (P s ) and (C p ) gives: 


T = p A d (V + w) (V + 2w) 

PC =T(V+w)=Tw 
s p o o 

solving for (w/V) gives: 

* . 0.5 /l + 4(y(i) 2 -0.5 


(7a) 

(7b) 


( 8 ) 
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or the thrust loss with forward speed is: 



( 9 ) 


During take-off and landing at velocities below 3w 0 , [note (w Q ) is computed from 
eq. 6a], the thrust loss can be approximated by: 


V 


When the CC blowing power (Pfc>) is extracted from the propeller shaft by an 
auxiliary compressor at an overall efficiency (vb)> then the remaining propeller 
power is reduced to (P 8 b) and the static induced velocity to (w Q b). 


| = 1 - 0.318 ( 


,0.678 


( 10 ) 


P, = 0.5 m . V. 2 (11) 

b J J 

P sb = P s - V \ < 12 > 

From the dependancy of (w D ) and (T 0 ) on the available power, the relative change 
in their magnitude with blowing is found from: 


w L 

ob 


w 

o 





(13) 

(14) 


When the propeller exhaust is vectored, the intake momentum (D m ) 
drag and only the outlet momentum is included in the thrust (Tb) 


is added to the 
or: 


D n c 

m _ 0.5 

Tob ~ w D b 


(15a) 



1 - 0.318 




0.678 


(15b) 


The blowing coefficient and CC lift augmentation reduce with increasing compressor 
pressure ratio and associated velocity (Vj): 
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(16) 
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If flaps are used for speed control then with a constant speed propeller, the 
compressor rpm, coupled to the propeller shaft, can also remain constant even 
during landing. This method is lighter and more reliable than using an auxiliary 
turbo-compressor to obtain constant CC blowing momentum. 


For jet aircraft, compressor bleed from either the last, or from an intermediate 
stage, produces a greater thrust loss than with shaft power extraction in a 
propeller aircraft. This is due to the loss of mass flow rate through the turbine 
and nozzle. An additional problem is that jet engine bleed air pressure decreases 
with RPM or throttle setting. To use this source for CC blowing air, one should be 
able to maintain a high power setting during approach to landing, which necessitates 
thrust vectoring. The associated engine thrust loss can be computed from cycle 
analysis on the T-S diagram. The magnitude of the component efficiencies and the 
temperature ratio’s for the compressor and turbine, determine the thrust loss 
associated with bleed air and forward speed. For example assume a compressor and 
turbine efficiency of 80% and 87.5% respectively, a turbine inlet temperature equal 
five times the ambient temperature and the compressor outlet temperature equals the 
square root of five times the ambient temperature. This results in a ratio of 
exhaust velocity to flight velocity given as a function of flight Mach number (M). 


Y 

Y = 2.37 M (17) 

The intake density and engine mass flow rate increase by ram compression to about: 

t-tJ r = 1 + 0.5 M 2 (18) 

m (V = o) 

resulting in an exhaust momentum increase with flight Mach number of about: 

Y = 1 + 0.5 M 2 (19) 

o 


The reduced static thrust (T 0 b) associated with bleed air, can be computed as a 
function of the bleed mass ratio (b = mj/m), note: limited to less than 20% for most 
engines. For example the thrust loss for high pressure bleed with 11 atm at the 
compressor final stage, is found to be 

T 

if*- = (1 - b) 2 “ (1 - 2b) (20a) 

o 

With 4 atm at the intermediate compressor stage, the thrust loss is found to be 

]ob = (1 - 0.5 b) (1 - b) « (1 - 1.5b) (20b) 

T 

o 

When thrust vectoring is employed the inlet momentum (m V = D m ) is added to the 
drag, and is about: 
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= 0.42 M 


(21a) 



The exhaust momentum variation with flight Mach number and bleed air is: 


T 

IjT = (1 + 0.5 M 2 ) (1 - b) 
ob 


(21b) 


The CC blowing momentum (Tj) obtained with the bleed air, depends on the 
temperature and pressure losses in the bleed air flow control valve and the ducting. 
Even when minimal losses are assumed, this results in a ratio of blowing thrust 
recovered to engine thrust lost of not more than: 

T . 

— r = 0.4 at 11 atm and 0.6 at 4 atm (22) 

o ob' 

This shows that extracting bleed air at the intermediate compressor stage is 
more efficient than at the final stage. The cycle analysis can be extended to 
include the effect of by-pass ratio when a turbofan is used. Experimental jet 
engine thrust loss data, including the effect of reduced bleed pressure at part 
throttle, are given by Hemmerly (1977). 


Ground Run in Landing and Take-Off 


The landing and take-off analysis presented here is based on the treatment by 
Kohlman (1981). During approach to landing, the descent angle is steep and the 
required forward component of the thrust is very small. To provide speed control, 
mechanical drag producing devices are needed unless the thrust can be vectored, 
for example using Pegasus type nozzles. In equilibrium flight, the minimum decsent 
velocity can be computed as a function of the available blowing momentum (Tj), from 
the requirement that the sum of the forces normal and along the direction of flight 
are zero. For CTOL aircraft with an approach speed equal 1.3 V s t,ali the approach 
angle of attack is found from 


“a = (23) 

u.3r 

For a CC powered STOL aircraft, the same safe approach angle of attack may be 
specified. If the flap size, deflection angle and effect on Cp 0 are specified then the 
dynamic pressure (q a ) on approach is found by iteration such that the descent 
angle (y a ), isolated from the two equilibrium force equations below, are equal. 
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y = cos 
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lii + 
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sin (e T + « q ) 


(25) 
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The resulting descent angle shall not be too steep, otherwise the pilot's 
visability is impaired or the vertical component of the approach velocity may exceed 
1000 ft/min. If such a problem occurs, then the flap or thrust vector angle must be 
reduced. An advantage of circulation control over other boundary layer control 
techniques is that high values of the lift coefficient can be obtained at moderate 
angles of attack. During a steep descent at an angle (y«) greater than the angle of 
attack, the aircraft attitude will be nose down, providing good pilot visibility. To 
obtain a flare-out before touch down one must be able to generate extra lift, by 
either having extra angle of attack or blowing pressure available. If one third of 
the lift is provided by each: the angle of attack, the CC lift augmentation, and the 

vectored thrust, then increasing the angle of attack for flare is only one third as 
effective as it is in a CTOL aircraft! During the optimum performance approach, 
there is no extra lift available, therefore the flare distance contribution to the 
landing may be ignored. The FAA specifies that after touch-down the deceleration 
rate is limited for passenger comfort to dV/dt - -0.5 g. If the deceleration rate is 
constant then the ground roll (s q ) is directly proportional to the minimum approach 
dynamic pressure or 


s (landing) = — ; 

p 1 dtJ 


(26) 


For take-off, the above equation shows that the ground run is likely just as 
sensitive to the minimum lift-off dynamic pressure (qiof) as ft i® ^h® acceleration. 
To minimize the lift-off speed one needs high blowing rates and high thrust angles. 
However, to maximize acceleration one needs all the obtainable thrust in the 
horizontal direction. Consequently optimum performance is obtained by delaying CC 
blowing and thrust vectoring to the moment of lift-off! Such a last minute 
configuration change increases the pilot's work load and reduces safety. Assuming 
there are no last minute configuration changes permitted, one finds that the optimum 
bleed air flow rate and thrust vector angle, become very important parameters in 
the take-off ground run distance optimization, more so than in the climb distance 
optimization. During the groundrun, on a low friction, level runway, the ground 
effect due to CC blowing may be ignored. At the start when q=o find C^=», and the 
lift increase due to CC should be ignored till C M reduces below 0.3. On the ground, 
the angle of attack (a g ) is constant prior to rotation and the ground distance is 
found by integrating the equation for the horizontal acceleration from (q w ) to (qiof)* 
Where (qw) i® the dynamic pressure of the head wind component and (qiof) i® the 
lift-off dynamic pressure. After rotation to ( a lof ) > the sum of the lift and vertical 
component of the thrust vector equals the weight. The take-off distance (sg) is: 


l lof 


s . = 


g 


1 

P g 


ti - /V ) dq 


(27) 
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'T fT 1 D 

^Jcos(e T +« g ) -f(l-l^rJsin(e T +« g ) -qC L /(")) - ^Vls) " ~W 


Near the end of the ground run, the rotation velocity is reached, and the pilot 
uses the elevator to increase the angle of attack by at least three degrees per 
second. Ideally, lift-off is achieved just prior to complete rotation to the maximum 
safe climb-out angle of attack. This ensures definite lift-off with an upward 
acceleration. To maintain a safe margin below the stall angle of attack, CTOL 
aircraft are specified to lift off at 20% above the stall speed or at an angle of 
attack: «i Q f = « stall/U* 2 ) 2 = 70% of the sta11 an « le * If the Bame ( a lof) angle of 
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attack is used on STOL aircraft, with powered lift and thrust vectoring, then the 
lift-off speed safety margin is less than 20% above the stall speed. When the 
lift-off angle of attack is specified, the lift-off dynamic pressure can be found from: 


0 = 



sm 


(6 T + “loP 


lof 

(W/S) 


aC T 

a, _ L 
lof — - + 
da 



+ C B 



(28) 


This is a quadratic equation in (qiof)^*® or 


'lof 



Maximum Climb Performance 

During the climb period needed to clear an obstacle, the thrust vector angle and 
the flap angle are not altered for safety reasons. The pilot only modulates the 
aircraft attitude and thus angle of attack to maximize the rate by which the 
maximum climb angle is reached. Because of the limited effectiveness of the angle of 
attack, it should as high as practical, while allowing an assymptotic approach to the 
maximum climb angle (r x ) at (« x ). A near optimum performance will be obtained 
when the angle of attack is reduced in proportion to the available acceleration, but 
never lower than (« x ). 


(« 


a = « x + 


lof 


\ / dV \ 

V ( dt> 


dV 

^~(at lift-off) 


(30) 


If the climb angle increases sufficiently fast 9 then the velocity does not 
overshoot the maximum climb angle value (V x ). In case it does exceed (V x ) the 
acceleration dV/dt will eventually become negative at which time (a x ) must be held 
constant. To apply this equation one must first compute the maximum climb angle 
parameters. 


From energy considerations, one may determine the steady state climb angle by 
equating the excess thrust power to the rate of increase in potential energy or 

[T b cos (e T + a) - D -D ] V = W.RC (31) 

The maximum climb angle is reached when: sin(y) = is maximum or 

sin = E(vrl c ° s <»t + »> - («1 - (irJha*. <32) 

When the engine thrust components are independent of flight speed, then the 
maximum climb angle occurs when the drag is minimum or (Cd q = Cpi). Without 
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powered lift (Cl) is proportional to («) and the induced drag coefficient is related 
to the lift coefficient squared. Then (Clx) an< * ( a x) are found from 


C, = -JcT. it e All = JcZ w e AR 
Lx Di Do 


(33) 


a 

x 


* c . 


"Lx - f 


'6 f 


*C, 


da 


(34) 


The maximum climb angle dynamic pressure (q x ) is obtained when the sum of the 
forces normal to V are zero or 

,2 

rfT i r 

2 , 


! <V r [(irl sin ( *T + “x J * S< C Lx / (I) 

Eliminate (r x ) from Eq. 32 and 35 by using: 


(35) 


cos 2 (7 ) = 1 - sin 2 
x 


= 1 - [iT sln (, t + “x> ‘ 20 DoV (I) - ^ 


(36) 


This gives a quadratic equation from which (qx) can be determined. In general, 
the thrust is not constant and the maximum climb angle values are found by 
iteration. For steady state climb, the sum of the forces parallel and normal to the 
direction of flight equal zero, they are respectively: 

o = [^] cos (« T + «) - ® ^ ~ sin (r) (37) 

o = (^) sin (« T + a) + - - cos (r) (38) 


At any flight speed, the tangent of the climb angle is obtained from their ratio. 


tan (y) = 


( T bl 
Iw J 

, ^ D 

cos (e T + a) - - 

- 

D 

m 

W 


1 

i 

* . 

sin (« T + «) + 

T 

w. 



(39) 


The maximum climb angle (r x ) is a function of angle of attack and found by 
iteration. Starting with V = Vi Q f and « = «iof» decrease (a) until tan(y) reaches a 
peak value, which represents a steady state climb. Higher climb angles can be 
found at higher velocities. Increment (V) and reduce («) by a small amount until 
the next steady state value for tan (7) is reached. This process is continued until 
the steady state value of tan (7) reaches a maximum, at which time the desired 
magnitudes of: V x , « x and 7 X have been obtained. The unsteady climb performance 
computation can now be performed. 
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Starting at lift-off with V=Vi G f and 7 = 0 , integrate the equations of motion until 
an obstacle of specified height has been cleared. The corresponding climb distance 
(s c ) must be added to the ground run to determine the total take-off distance. With 
powered lift, all the aerodynamic coefficients such as: lift, blowing momentum and 
drag are functions of the dynamic pressure (q). The second most important 
parameter is the climb angle ( 7 ) or time (t). In general only time is a montonic 
increasing variable and both acceleration component equations (# and x to V) must 
be incorporated in the numerical integration. 


dt " g [(w*) COS (0 T + a) " W " "W ” Sin(r) . 
dt 2 = g [(ir) Sin (e T + a) + W " cos(7) . 


(40) 

(41) 


The initial value of the derivative (dy/dt) is zero if the climb-out is started 
with a = «i 0 f. However in practice the lift-off rotation may be continued to start 
climbing with («) slightly greater than («l 0 f) and (dy/dt) finite. If ( 7 ) increases 
monotonically, then the ratio's of the two acceleration equations may be used to 
eliminate the variable (dt) the result being a single equation in terms 
of (dy/dq): 


dz _ . 

[ T b) 

lw J 

sin (e T + a) + jj - cos ( 7 ) 

dq_ 2q 

f T b] 
lw J 

n D 

cos (e T + «) - - ® - sin(y) 


(42) 


This technique is not applicable to a ski-jump take-off where ( 7 ) first 
decreases before it increases. Then one must integrate Eqs. 40 and 41 with respect 
to time. Once the values of the (V), ( 7 ) are found at each time increment, the 
corresponding horizontal and vertical distance increments are found from: 

dx = dt [V cos ( 7 )] (average) (43) 

dy = dt [ V sin ( 7 )] (average) ( 44 ) 

The integration can be terminated as soon as the vertical distance (y) exceeds 
the obstacle height (h). The corresponding climb-distance is found by interpolation: 


S c = x “($) ( y “ h > 
and the total take-off distance is 


(45) 


s = s + s 
C g 


(46) 


Because the dynamic pressure (q) is inversely proportional to the wing loading 
(W/S), all the distance dimensions can be computed as the ratio s/(W/S), when the 
the obstacle height is specified in terms of h/(W/S). The velocities are then 
found in terms of V/(w/s)0>5. 
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For CTOL air craft the maximum climb angle and the difference between V x and 
Vjof are small. In such a case the steady state climb is reached quickly and 
prior to clearing the obstacle. For STOL aircraft, the maximum climb angle as 
well as the difference between V x and Vi 0 f are large as shown in fig. 10. In 
that case, the obstacle may be cleared prior to reaching the steady state 
conditions. Due to the limited rate of increase in the climb angle it is 
possible to overshoot (V x ) prior to reaching the steady state climb angle (r x ). 
When the pilot modulates the angle of attack as suggested in eq. 30, then the 
time required to reach the steady state climb angle (r x ) is minimal. Because 
(r x ) decreases with bleed air, it is advisable to decrease the blowing rate 
gradually during climb-out. 


Example of a Take-Off Performance 

To determine the minimum runway length for a CC STOL aircraft, the 
designer must compare the performance over a wide range of bleed air mass 
flow ratio's and thrust vector angles. The associated lift and blowing 
coefficients vary drastically and should not be used as design input variables. 
Suitable input parameters are: take-off thrust to weight ratio, thrust vector 

angle, blowing momentum as a function of bleed air mass ratio, thrust loss as 
a function of blowing momentum, and lift proportionality constant (Cb). All 
other input parameters are similar to those needed for a CTOL aircraft. A 
sample calculation was performed for an aircraft with the following 
characteristics: 

Tq/W = 0.5, (e AR) = 5, aCi/ a « = 4.74, Cfi = 5.6 «g = 3‘, cq 0 f = 11*> 6 f = °* 

blowing air supplied at 11 atm, and duct momentum loss 20%. The distances for 
both the ground roll and the climb distance are given as a function of the 
wing loading (W/S). This requires the obstacle height to be specified in the 
same units. Here (h = 0.5), this means (h = 35 ft) for a wing loading of 70 
psf. The ground run distance is shown in fig. 11 to reach a minimum at 5.5% 
bleed air and a thrust vector angle of 24 degrees. The minimum ground 
distance in ft is 15.3 times the wing loading in psf. This is less than half the 
distance required without CC or thrust vectoring. At this low thrust to 
weight ratio the effect of thrust vectoring by itself is about 20%. The loss of 
thrust and acceleration associated with high bleed ratios and thrust vectoring 
are clearly noticeable. In fig. 12, both the climb distance to clear an obstacle 
and the total take-off distance are shown. The penalty of not using the 
optimum bleed air ratio and thrust vector angle is not as critical for the climb 
distance as for the ground distance. This would be different if the bleed air 
extraction were delayed till the moment of lift-off. The usual obstacle height 
is 15 meter or 50 ft. If those distances were used, the climb distance would 
be closer to the ground distance. The trimmed lift coefficient, shown in fig. 
13 varies widely with blowing rate from below 1.0 to above 4.0 and for the 
minimum take-off distance, it is significantly below the maximum obtainable 
value. The corresponding blowing coefficient varies with the lift-off dynamic 
pressure as shown in fig. 14. Like the lift coefficient, it varies over a wide 
range of values. Practical values are limited to below 0.3, above that value 
the lift augmentation ratio equal A CL/(Cp)^*5 becomes less effective. 
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Conclusions 


Designing with CC requires a good insight in CC airfoil aerodynamics, the 
relation between engine thrust loss as a function of blowing air generation 
and the lift loss associated with trimming the pitching moments. The large 
number of design variables involved, necessitates numerous performance 
calculations were each parameter is varied over a wide range. The linear 
equations presented in this paper facilitate such an analysis. The thrust loss 
equations should be modified in accordance with the characteristics of the 
engine under consideration. Because cost, weight, simplicity and reliability 
are over-riding factors in the optimization process the minimum achievable 
take-off and landing distances may not represent the optimum design. The 
finally selected percent bleed air and duct size are probably lower than those 
corresponding to the minimum runway length. The partial derivatives of these 
parameters, such as (ds/db), should be plotted as a function of the high lift 
system cost or weight, to arrive at an optimum design. 

In a comparitive study with other high lift techniques one must realize 
that even though the power efficiency of CC is very attractive, the magnitude 
of the locally obtainable lift coefficient, is less than six. Desirable flying 

aspects are good pilot visability and effective DLC. Undesirable 

characteristics are sudden wing stall, and required in flight stowing of the 
rounded trailing edge. From the WVU CC flight test, it was found that drag, 
with the forward folding CC flap in the downward position, was so large that 
at full throttle, the sink rate was 1000 ft/min. Because the flap folding 
maneuver could be completed within four seconds, the change in speed and 
pitch attitude was acceptable. 

An improvement in the stowable flap design was recently patented, see 
figure 15. Here an independently operating Fowler flap is equipped with a CC 
rounded trailing edge, when needed. Safety and descent control are enhanced 
by the availability of conventional flaps in conjunction with CC. A separate 
air duct, which doubles as a rounded Coanda surface, can be swung in 
position, like an agricultural airplane spray bar. The available internal 
ejector suction provides the attachment force between the flap trailing edge 
and the CC air duct. The duct pivot point can be positioned anywhere to 
accommodate duct stowing at an optimum location in the wing. The drag 
associated with stowing only the blowing duct, is much lower than with a 
forward folding flap. 


583 



References 


1. Boasson, M., Loth, J.L., "Combination of Suction and Tangential Blowing in 
Boundary Layer Control" 27th Israel Conf. on Aviation and Astronautics, 
Feb. 1985. 

2. Carr, J.E. "Blended Blown Flaps and Vectored Trust for Low Speed Flight", 
J. Aircraft, Vol. 23, No. 1 1986 p. 26-31. 

3. Davidson, I.M., "Aerofoil Boundary-Layer Control System", British Patent 
No. 913754, 1960. 

4. Englar, R.J., "Circulation Control for High Lift and Drag Generation on 
STOL Aircraft. Journal of Aircraft. Vol. 12, No. 5, May 1975, pp. 457-463. 

5. Englar, R.J., Trobaugh, L.A. and Hemmerly, R.A., "STOL Potential of the 
Circulation Control Wing for High-Performance Aircraft", J. of Aircraft, Vol. 
15, No. 3, March 1978, pp. 175-181. 

6. Englar, R.J., "Low-Speed Aerodynamic Characteristics of a Small, Fixed 
Trailing Edge Circulation Control Wing Configuration Fitted to a 
Supercritical Airfoil. DTNSRDC/ASED-81/08, March 1981. 

7. Hemmerly, R.A., "An Investigation of the Performance of a J52-P-8A Engine 
Operating Under the Influence of High Bleed Flow Extraction Rates", DT 
NSRDC ASED-387, August 1977. 

8. Kind, R.J. and Maull, D.J., "An Experimental Investigation of a Low-Speed 
Circulation Controlled Aerofoil", Aeronautical Quarterly, Vol. XIX, May 10, 
1968, pp. 170-182. 

9. Kohlman, D.L. "VSTOL Airplanes", Iowa State U. Press, 1981. 

10. Loth, J.L., "Some Aspects of STOL Aircraft Aerodynamics", SAE Paper No. 
730328, Business Aircraft Meeting, April 1973. 

11. Loth, J.L., Fanucci, J.B., Roberts, S.C., "Flight Performance of a Circulation 
Controlled STOL Aircraft", Journal of Aircraft, Vol. 13, No. 3, March 1976, 
pp. 169-173, also AIAA Paper 74-994, August 1974. 

12. Loth, J.L., Boasson, M., "Circulation Controlled STOL Wing Optimization," J. 
Aircraft, Vol. 21, No. 2, 1984, pp 128-134. 

13. Roberts, S.C., "WVU Circulation Controlled STOL Aircraft Flight Tests", 
WVU Aerospace Technical Report No. 42, July 1974. 

14. Williams, R.M. and Howe, H.J., "Two-Dimensional Subsonic Wind Tunnel Tests 
on a 20 Percent Thick, 5 Percent Combined Circulation Control Airfoil. 
NSRDC TECH. Note AL-176 (AD877-764) 1970. 


58M 



ORIGINAL PAGE TS 
OF POOR QUALITY 


STOL CHORD 02=1.2 Ci 
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LE DROOP 


3“ DECOMPRESSED 
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1 kG /SEC AT 28 PSIG 


2 kG /SEC 


Figure 1 - CC airfoil used on the 
WVU CC Demonstrator STOL 
aircraft. 
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Figure 2 - Hollow flap of the WVU 
CC wing with supersonic jet 
ejector. 



Figure 3 - WVU STOL aircraft with 
aerodynamic and structural fences 
needed to transfer loads past the 
CC duct stowing cavity. 


Figure 4 - CC air supply with DLC 
sliding air dump valve and aileron 
linked splitter vane for roll control 
by differential blowing. 


CL/R.VE FIT • C^_ =4.74 («<-*- 0.02.8 aapm/v) + 6i 

W/*0 

Figure 5 - Empirical curve fit with 
Eq. 2| applied to wind tunnel data 
on A-6 model with tail-off. 
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Figure 7 - Linear CC airfoil lift 
increase versus velocity ratio. 


Figure 6 - Non-linear CC airfoil lift 
versus blowing pressure. 
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Figure 8 - Effect of velocity ratio 
on lift increase at constant power. 
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Figure 9 - Example of reduction in 
blowing momentum due to duct 
loss. 
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Figure. 10 - Reduction in lift-off 
and best climb angle speed with 
CC. 
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NOMENCLATURE 


Ad 

AR 

b 

C B 

CC 

C D 

Cl 

Cp 

cm 

c 

Djn 

e 

f 

i 

h 

M 

m 

Ps 

q 

RC 

s 

s 

t 

T 

V 

W 

w 


Propeller disk area 
Wing aspect ratio 

Ratio of bleed air to jet engine mass flow 
Lift proportionality constant for CC 
Circulation Controlled abbreviation 
Total airplane drag coefficient 
Lift coefficient 

Coefficient of preformance for a propeller 

CC blowing coefficient 

Wing chord in cruise mode 

Engine inlet air flow momentum 

Span-wise loading efficiency 

Runway friction coefficient 

Acceleration of gravity 

Take-off obstacle height to clear or CC blowing slot height 
Free stream Mach number 
Mass flow rate 
Shaft power 

Free stream dynamic pressure 

Rate of Climb 

Wing area in cruise 

Horizontial take-off or landing distance 
Temperature or time 

Thrust, except with thrust vectoring is only exhaust momentum 
Velocity, without subscript means flight speed 
Aircraft weight 

Velocity increase induced at the propeller disk 


Subscripts 


a 

b 

c 

e 

g 

i 

j 

lof 

o 

T 

w 

X 

Greek 


a 

7 

6 f 

V 

e 

p 


approach related parameter 

parameter with CC blowing 

climb related parameter 

exhaust parameter 

ground distance 

induced drag 

blowing jet parameter 

parameter related to lift-off speed 

parasite drag or reference static thrust 

thrust related parameter 

head wind component 

maximum angle of climb parameter 


angle of attack 
flight path angle 
flap angle 
efficiency 

thrust vector angle or flap hinge locater angle 
air density 
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